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^sj ■ A compact gas cloud G2 1 is predicted to reach the pericentre of its orbit around the super 
massive black hole (SMBH) of our galaxy, Sagittarius (Sgr A*), by summer 2013. This event 
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will give us a rare opportunity to observe the interaction between SMBH and gas around it. 

rS ! 
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We report the result of the fully three-dimensional simulation of the evolution of G2 during 
the first pericentre passage. The strong tidal force by the SMBH stretches the cloud along 
its orbit, and compress it strongly in the vertical direction, resulting in the heating up and 
flaring up of the cloud. The bolometric luminosity will reach the maximum of 100 L & by July 
2013. This flare should be easily observed in the near infrared. 
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At present, the activity of Sgr A* seems to be in the low-state, with the X-ray luminosity of 
10 33 erg s -1 ^. There are evidences of past activities^, where the luminosity had reached as high 
as ~ 10 40 erg s^ 1 . Recently, rapid flaring from Sgr A* was observed in various 

bandsPHS 1 . xhus 

it is quite important to understand how these large variations in the luminosity took place. One 
possibility is the intermittent supply of gas in the form of high-density clouds. 

The compact cloud G2 1 might be offering us an unprecedented opportunity to study such 
an interaction of a gas cloud and SMBH. While the formation mechanism of the cloud is under 
debateQIEE103], we k n0 w that the orbit of G2, determined from observations starting at 2004, is highly 
eccentric, and G2 will reach the pericentre of its orbit in the summer 2013, with the pericentre 
distance of only 270 au. 

Schartmann et alJ^ studied the evolution of G2 using high-resolution AMR calculation in 
two dimensions. In their calculation, the cloud loses the kinetic energy during the pericentre pas- 
sage due to the ram pressure from the hot atmosphere around the SMBH, and gas accretion to 
SMBH starts in early next year, and continues for coming several decades. 

However, it is not clear if the two-dimensional calculation is appropriate or not. The cloud 
should experience strong compression in the direction perpendicular to the orbital plane, due to 
the tidal force from the SMBH. If we neglect the effect of the pressure, the cloud would become 
razor-thin near the pericentre. Therefore, the ram pressure might be much less effective as the drag 
force to the cloud. Moreover, the compression heats up the cloud, which should emit radiation. 
Thus, there might be an increase in the luminosity of the cloud, solely due to this compression. 
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In order to study these effects, a fully three-dimensional simulation is necessary. We performed 
such three-dimensional calculations, in which the compressional heating and radiative cooling of 
the cloud are consistently taken into account. 

Figure 1 shows the evolution of the three-dimensional structure of the cloud from A.D. 2006 
to A.D. 2014. The initial orbital elements of G2 were determined from the observation^. We 
adopted the compact cloud scenaricP-R For this model, we used the atmosphere model which is 
similar to that used in the previous 2D calculation™!. At A.D. 2006, the simulated G2 was nearly 
spherical. By 2012.2, it is stretched in the orbital plane and compressed in the vertical direction. 
This stretch of G2 has already been observed 1 . When the cloud passes the pericentre, its thickness 
reaches the minimum. Due to this strong vertical compression, the gas density increases by more 
than two orders of magnitude. 

In the last two panels (panels e and f), we can see the "bridge" between the central SMBH 
and the head of the cloud. This bridge indicates that there is a flow of gas from the head of the 
cloud to the SMBH. However, the amount of the gas in this bridge, and the resulting accretion rate 
to SMBH, are very small. 

Figure 2 shows the time evolution of accretion rate. We show the results from two runs out 
of the three runs which we have tried. For the first run, we used the cloud model the same as that 
used in Burkert et al^ and Schartmann et al. except that we solve the dynamical evolution of 
the hot ambient. Whether or not such high-density, high-temperature atmosphere actually exists 
near SMBH is an open question. In order to test the importance of the assumption on the hot 
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Figure 1 : Time evolution of the structure of G2 in the three dimensional simulation for the period of A.D. 
2006.26 to A.D. 2013.82, tagged with a to f and a' to f. The density distributions in the xy planes (a to f 
) and the selected planes (a' to f) are shown. Here, the xy plane is corresponding to the orbital plane of 
G2. The thin solid lines in the panels for xy planes show the location of the plane for which the density 
distributions are shown in the panels in the corresponding panels for selected planes. The top two panels 
shows the distribution at A.D. 2013.46. The middle and bottom rows show the view in xy and selected 
planes, respectively. The panel g depicts part of the orbit of G2 (solid curve) and the positions of the xy 
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atmosphere, we performed two additional simulations. In the second run, we reduced the gas 
density by a factor of 10 (red dashed curve in the figure), and in the third run we eliminated the 
atmosphere altogether. 

In the case of the standard run, the accretion rate of the gas from the cloud reaches the peak 
value of ~ 1O -7 M year -1 at A.D. 2014, and then decreases exponentially. The accretion rate at 
A.D. 2030 is one order of magnitude smaller than the peak value. 

These behaviours are quite different from those in two-dimensional calculations 1 , in which 
the accretion rate is nearly constant due to the strong ram pressure. In our three-dimensional simu- 
lations, the height of the cloud is reduced to 1/100 of the initial value at the pericentre, resulting in 
the decrease of the effect of the ram pressure by a similar factor. Even if we assumed a relatively 
high density for the halo gas, it is probably difficult to observe the change of the activity of Sgr A* 
due to this additional accretion, since it is one or two orders of magnitude lower than the estimated 
accretion rate of the hot ambient gas, which is 10~ 6 ~ 10~ 5 M year -1 . 

The overall evolution of the accretion rate depends strongly on the assumed density of the hot 
atmosphere. When we reduced the hot gas density by a factor of 10, the accretion rate decreased 
by the same factor. In the run with no halo gas, no gas is accreted to the SMBH. 

Figure 3 shows the evolution of the total luminosity of the cloud G2, for three different 
models for the hot atmosphere around the SMBH. We neglect the emission from the gas absorbed 
by Sgr A*, because the accretion around Sgr A* is expected as radiatively inefficient accretion flow 
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Figure 2: Time evolution of the accretion rate of the gas to the SMBH. Black and red curves are the 
results of simulations with the standard high-density hot atmosphere (Runl) and with low-density 
atmosphere (Run2), respectively. Thick curves show accretion rates of the cloud gas, and thin 
curves show those of the halo gas. The vertical line denotes the epoch of 01/09/2012. We used the 
accretion radius of the SMBH of 30 au, which is 400 times larger than the real horizon scale of 
SMBH, 0.085au, and is 10 times smaller than the pericentre distance. When we excluded the hot 
ambient, there was no accretion to the SMBH. 
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Figure 3: Total luminosity as a function of time. Panel a shows the time evolution of the luminostiy 
from A.D. 1995 to A.D. 2030, while panel b shows it from A.D. 2011.5 to 2014.5. Solid black, 
red dashed and green dot-dashed curves are the results from simulations with the standard high- 
density hot atmosphere (Runl), with low-density atmosphere (Run2), and no atmosphere (Run3), 
respectively. 

(RIAF). The peak luminosity and its duration is practically independent of the assumption for the 
hot atmosphere. This result is quite natural, since the peak luminosity comes mainly from the tidal 
compressional heating of the cloud and has nothing to do with the interaction with the atmosphere. 
On the other hand, the interaction with hot atmosphere keeps the luminosity high for years before 
and after the pericentre passage. 

We can estimate the total amount of the energy generated by the tidal compressional heating 
in the following way. For simplicity, let's assume that gas moves freely until it reaches the equato- 
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rial plane, where it converts all the kinematic energy of its vertical motion to thermal energy. The 
vertical velocity is given by V v = V a tan i, where V a is the velocity at the ascending node and % is 
the inclination. If we assume that the cloud was still spherical around year 2000, we have i ~ 5° 
and V v ~ V a tanz = 440 km s -1 for a gas element 125 au from the orbital plane at the year 2000. 
By integrating the energy over the spherical cloud of radius R, we have 

dt V 125au / V 3M ffi/ V r J 

* 1.1 xW % Vf— V^W- 1 . (1) 

V5 degree y> \3Mj \ r ) 
Here, R is the radius of the cloud when the cloud is on the minor axis of its orbit. We assume 
that vertical velocity is zero at that moment. In other words, we assume that the ascending node 
coincides with the pericentre. The inclination i is that of the gas element with the largest distance 
from the equatorial plane. The duration time of energy release r is about one year, as we can see 
from figure 3. 

This simple estimate is in good agreement with the total amount of radiation in our detailed 
simulations. On the other hand, the total amount of radiation depends on the height and vertical 
velocity structure of the cloud. For example, if the compression velocity is zero at the apocentre, 
the inclination would be smaller by a factor of a few, resulting in the decrease of the total luminosity 
by about one order of magnitude. 

Figure 4 shows the distributions of temperature and brightness, for selected moments. The 
main body of the gas cloud is heated up by the compression, but the temperature remains near 
10 4 K. This is due to very efficient radiative cooling of ionized hydrogen through recombination 
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Figure 4: Time evolution of the gas temperature and integrated surface luminosity from gas with a 
temperature range of 8000 - 15000 K from A.D. 2006.26 to A.D. 2013.82, tagged with a to f and 
g to 1. The face-on view (the xy plane) is assumed and each panel shows a 900 au x 900 au region. 



lines. The most luminous region is several hundred aus in size. They are most likely to be observ- 
able from the outside in the near infrared, which does not suffer from the dust extinction effects 
significantly. 
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We have performed fully three-dimensional simulations of the evolution of the G2 cloud. 
Our result differs from the result of previous two-dimensional simulations in (a) strong vertical 
compression leads to the heating up and flaring up of the cloud at the first pericentre passage, 
and (b) because of this compression, ram pressure drag from the hot atmosphere is ineffective in 
removing the energy and angular momentum of the cloud. 

In our standard model, the peak luminosity would reach 100 times the solar luminosity. 
The luminosity depends on the assumed internal velocity structure of the cloud, and thus might 
be fainter by one order of magnitude. Since the peak luminosity is from the tidal compressional 
heating, it does not depend on the assumption on the structure of the hot atmosphere around the 
SMBH. We therefore believe that our prediction is pretty robust. 

We expect that the increase of the luminosity of the cloud would have already started and 
should be observable in near infrared now. In addition, the increase of the vertical velocity is 
probably observable as the line broadening. Since the vertical velocity should strongly depend on 
the position of the gas on the orbit, it is very important to measure the variation of the velocity 
profile both in space and time. 

Detailed comparison between high-accuracy three-dimensional calculations and observation 
will help us to understand the nature of the cloud and how it will interact with the SMBH. 
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